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Are sterile neutrinos consistent with clusters, the CMB 
and MOND? 
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INTRODUCTION 



MilRrom's Modified Newtonian Dynam i cs (MOND, see 
lMilErom"l984 ISanders fc M cGaugh' 2002^: lBekenstein|[200^ : 



iMilgrpin 20081 ) as the weak acceleration limit of Einstein's 
' general relativity is arguably consistent with a wide range 
' of evidence from astronomical systems from the orbits of 

the planets in the solar system (where Newto n's gravity suf- 

■ fices withou t dark matter) to dwarf galaxies ( Milgrom 19951; 
\ lAngusl200i ') and globular clusters (Angus & McGa ugh 2008, 

■ in prep) of the Milky Wa y, tidal d warf galaxies (jMilgronj 
' l2007bl : iGentile et al] |2007| . see also 'Bournaud et al ' 



low surface brightness ga laxies (McGaugh & de Blok 19981 : 
iMilgrom fc Sander s'20071 and high s urface brightness galax- 
ies llFamaev fc Bhinev..2005. : .Sanders fc Noordermeer] 20071: 
McGaughT 20081 ) including giant ellipticals ( Angus et al.l 



ABSTRACT 

If a single sterile neutrino exists such that m^, lleV, it can serendipitously solve 
all outstanding issues of the Modified Newtonian Dynamics. With it one can explain 
the dark matter of galaxy clusters without influencing individual galaxies, match the 
angular power spectrum of the cosmic microwave background and potentially fit the 
matter power spectrum. This model is flat with H.^^ ~ 0.23 and the usual baryonic 
and dark energy components, thus the Universe has the same expansion history as 
the KCDM model and only differs at the galactic scale where the Modifled Dynamics 
outperforms ACDM signiflcantly. 

iGentile et"aLl |2007| . see also iBournaud et al]|2007^ . Never- 
theless, it is the generally accepted model. 

One might find this astounding except for the fact 
that it has resounding success at cosmological scales where 
MOND predictions are sketchy. Furthermore, clusters of 
galaxies r equire large quantities of dark mat t er in MOND 
(jSandcrs 2003: 'Pointocouteau & Silkl l2005l: IC lowc et al.1 
[2006; Angus et ah, 2007; Sanders 200 7i; lAngus et al. .2008]^ 
This might sound like a contradiction, but it is perfectly 
sensible as long as an obvious contraint is satisfied i.e. that 
any dark matter in MOND has a free streaming scale that 
is larger than galaxies and the density is low. Otherwise the 
dark matter would manifest itself in ordinary galaxies which 



2007I) . 



2008 



Not only are the dynamics of the galactic systems well 
matched from the MOND prediction , but they all fall pre- 
cisely on the Tu Uy-Fisher relation (|McGaugh ct al.| l2000l : 
lMcGaughir2005bl ) which correlates total enclosed mass and 
the fourth power of the asymptotic velocity, unless they are 
satellites of a larger galaxy. 

The only alternative theory uses cold dark matter 
(CDM) particles without any experimental motivation in 
massive, triaxial halos to provide the additional grav- 
ity needed to boost the rotation velocities of the sys- 
tems with an acceleration discrepancy. There are sev- 
eral well-documented and as yet unresolved issues for 
this framework at the scales of galaxies such as the fine 



McGaughll2005al). the cusp problem (de Blok fc McGaugh, 


19981: McGaugh fc de Blokl 19981: iGnedin & Zhad I2OO2I: 


Gentile et al.1 12004 iGilmore et al. 


20071), the missing 


satellites problem (IKlvDin et al.l 199!: 


: Moore et al. 


|l999|) 


and more recently tidal dwarf galaxies (MilgromI 


2007bl: 



would d estroy the consistency of the baryonic T uUy-Fisher 
relation (|McGaugh et aLlbOOrt 'McGaugh'2005bi). 

It was postulated (ISandera ,2003, 2007 ) that the active 
neutrinos (at = 2eV; very close to the experimental 
upper limit of 2.2eV) can provide the dark mass of clus- 
ters. Neutrinos conform to certain scaling relations in clus- 
ters such as the proportionality of the electron density in 
the cores of clusters to T^/^. Nevertheless, recent s tudies 
l|Pointecouteau fc Silkl l2005l : lAngus et al] |2007| . I2OO8I ) have 
shown that, even under very favourable circumstances, a 
second species of dark matter would be necessary to explain 
the dynamics of the central lOOkpc of clusters and groups 
of galaxies, which more or less rules them out the active 
neutrinos as good candidates. 

Of course, there is no limit on the dark matter being 
baryonic since the necessary dark matter in clusters of galax- 
ies is at most a few percent of the big bang nucleosynthesis 
(BBN) baryons, of whi ch only 20% or so ar e observed at 
low redshifts (Silk 2007: IMcGaugh et al.ll2007l '). the remain- 
der is pre s umed to exist in a warm-hot inter galactic medium 
iBregmanI l|2007l '). This led iMilgromI to propose the 

dark matter in clusters to be a cold, molecular gas of ~ 
Jupiter mass. These are naturally difficult to detect, but 
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might have the se rendipitous fortune of resolving the cool- 
ing flow problem (jFabian et al.lll994l ). 

Unfortunately, even if the cluster dark matter problem 
were resolved, there remains the issue of cosmological dark 
matter. Put simply, there is compelling evidence that the 
Universe consists of a form of dark energy (like a cosmologi- 
cal constant) that forc es the expansion of the Universe to ac- 
celera te at late times (jPerlmutter et al]ll99 9: Schm idt et alJ 
Il998|). However, we have no idea what this dark energy 
IS IPiaferid lioosh from a particle physics point of view, 



although perhaps the coincidence between Uo and c Hp or 
c(A/3)^/^ is a strong indication (Milgrom 2002. |2008| ). 

With the presence of this dark energy, in order for the 
Universe not to expand too rapidly, there needs to be some 
form of matter independent of the well fixed quantity of 
baryons to endow the Universe with additional inertia. This 
additional matter serves several purposes: it allows for large 
structures to form more rapidly out of the expanding Uni- 
verse which is shown by the matter power spectrum at large 
scales l|Tegmark et al.l 120041 ). Also, it drives the collapse of 
the photon-baryon fluid to form fluctuations in the cosmic 
micro wave background ( CMB) at well measured angular 
scales (IWhite et al.lll994l ) and it gives the correct distance- 
redshift relation (expansion history). 

The underlying theory of MOND is still unknown, also 
it is only a classical framework, so a huge effort was made 
to extend MOND to the relativistic reg ime. In 2004, a gi - 
ant leap was made in this direction by iBekensteinI (|2004l ) 
and others have taken to thrashing out the predictions of 
other MOND insp i red relativistic theories ([Sanders! l2005l : 

■ M- 



ISkordis et ahlbood : IZlosnik et al.ll20o'^ . l2007al lbD. Sadly, the 
predictions for cosmology are not clear and there seems to 
be too much freedom, in contradiction to the absolute pre- 
dictiveness of MOND in galaxies. 

For this reason, I show here the predictions of coupling 
MOND wit h sterile neutrino dark matter using the ansatz 
employed bv lMcGaugiil (|2004i ') when matching the CMB with 
MOND i.e. that no MOND effects are present before recom- 
bination. A simple argument supporting this is that at a 
redshift of z ~1080 the angular diameter distance to recom- 
bination Da = 14Gpc and the angular scale, 9, of the first 
(and largest) peak is 1° or 0.017rad. So the physical size of 
the first peak r = 6Da is ~240Mpc. Since the average over- 
density 5 is only 1 part in 10^ of the critical density pc{z), 
the typical gravities at a radius r from the centre of one of 
these overdensities is g — G5M{r)r~^ — ^GS pc(z)r where 

= ^iS- and H{zf = Hll^mil + zf ^Ua. Compiling 
all this gives 



g{r) ~ ]^5Hl + zf + Oa] r 

At r = 240Mpc, 



(1) 



It is often forgotten when looking at MOND cosmol- 
ogy that no cold dark matter exists in MOND. Therefore, 
we must relax many of the constraints that are set by 
CDM cosmology. The most important and obvious one is 
that there is now a large gap in the energy-density bud- 
get since CDM is not present and it is perfectly reason- 
able to fill this gap with hot dark matter like neutrinos. 
The constraints on neutrino masses, for which cosmology 
is still the most stringent, must be reanalysed in light of 
MOND. Still, the empirical evidence from su pernovae data 
(|Schmidt et al.lll998l : IPerlmutter et al.lll999l ') strongly sug- 
gest the universes expansion is accelerating owed to the ex- 
istence of dark energy, Qa- Furthermore, the baryon bud- 
get is strongly constrained by well understood physics to be 
around Qbh'^ 0.015 -0.025 (iBoesgaard fc Steigmanlll985l : 
iBurles et al.ll200ll : lMcGaughll2004h . but this still leaves a 
large amount of latitude in the energy budget for DM. 

Any DM, however, must be compatable with clusters 
of galaxies, the well understood lack of DM in galaxies in 
MOND and the anisotropies in the angular power spectrum 
of the CMB. The best candidates for such hot DM are neu- 
trinos. 



2 NEUTRINOS 

2.1 Active Neutrinos 

The three active neutrinos {i/fj,, i/e and i/t) from the stan- 
dard model of particle physics have been shown to mix be- 
tween flavours by atmospheric and solar n eutrino experi- 
ments (I Ahmad et al.ll200ll : lAshie et al.ll2004l '). However, the 
exact masses of the three active neutrinos are not yet known, 
only their squared mass differences. Nevertheless, the masses 
of all three are kno wn to be less than 2 .2eV from the Mainz- 
Troitz experiments l|Kraus et al.ll2005l ). 

The maximum density that a neutrino species can pro- 
duce after gravita tional collapse is given by the Tremain- 
Gunn limit (jTremaine fc Gunnlll979l ), 



9 ' 



i • 10"^ ■ 7.1 X 10"^ [0.27 ■ 1081^ + 0.73] ■ 240Mpc ~ 570ao, 



where ao — 3.6(kms~^) pc~^ is the MOND acceleration 
constant. Typical accelerations so many times greater than 
Oo are completely unaffected by MOND gravity and there- 
fore no MOND effects should influence the CMB. However, 
as z drops, so does pdz) and thus peculiar accelerations 
can slide into the MOND regime. Thus, the matter power 
spectrum can be affected by MOND. 



7 X 10-^Mqpc- 



\lkeVj 



\2ev) 



(2) 



for each of the three species. Thus, the density is greatly 
dependent on the mass of the neutrinos. However, groups 
and clusters of galaxies have dark matter that is much denser 
than can be produced by the active neutrino s even at the 
maximum mass of 2.2eV l| Angus et al.|[200^ . If the dark 
matter is i ndeed a neutrino li ke species, it must be heavier 
than 8eV (| Angus et al.|[200^ . There is a further problem 
with neutrinos at 2.2eV in that the contribution they make 
to the energy density of the Universe is given by 



f^i, = 0.0205m^, 



(3) 



meaning that at 2.2eV the three neutrinos make a 13.6% 
contribution to the energy density of the Universe, but the 
maximum density is relatively low (see Eq[2]). Such a huge 
contribution would be easily detectable in the angular power 
spectrum of the fluctuations in the CMB as shown for this 
example in Fig (2] Therefore, the active neutrinos are a very 
poorly motivated candidate. 
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2.2 Sterile Neutrinos and the CMB 

As mentioned above, the three active neutrinos are known 
to have mass. Another oddity arising from this is that the 
active neutrinos are solely left handedly chiral, whereas all 
other fermions are ambidextrous. The easiest way to incor- 
porate this into the standard model of particle physics is 
to introduce a right handed "sterile neutrino" . In addition, 
they are not simply aesthetically pleasing, the introduction 
of a single sterile neutrino was preferre d from analysis of the 



Miniboone experiment by Giunti fc Laveder ( 2007) (see also 
lAguilar et al.ll200ll : lMahoni fc SchwetjboOTl ') with a mass in 
the range 4eV-18eV to explain the disappearance of electron 
neutrinos from the beam at low energies. 

In the simplest model, if the mixing angle of the ster- 
ile neutrino is low enough, then thermalisation in the early 
Universe can balance the abundance of the sterile and active 
neutrinos. In this case, the cosmological density is exactly 
related to their mass, as for the active ones (EqO. 

With the hypothesis that all the DM in MOND comes 
from a single sterile neutrino, we used the free ly available 
CMB anisotropy code CAMB (|Lewis et al.ll2000l ) and incor- 
porated it into a minimisation routine c omparing with 



the d ata from the WMAP5 data release (iDunklev et al 



20081') and the ACBAR 2008 data release (|Reichardt et al 



20081 ). We allowed variation of fit, f^i^si ^a, dus/dlnk, r. Ho 
and fixed the Universe to be flat meaning Qa = 1 — Sib — ■ 
Obviously, in this MOND inspired model there is no 
CDM by definition, but since the CDM model works well 
at producing the CMB anisotropies, we began the search 
by simply transferring flcdm to fi^^ . Furthermore, the 3 ac- 
tive neutrinos are taken, for simplicity, to be massless. As 
discussed later, it is not feasible to have a pair of very mas- 
sive (> O.beV) sterile neutrinos because splitting the Q^^ 
between two or more neutrinos reduces the available mass 
to each neutrino thus detrimentally lowering its Tremaine- 
Gunn limit {p^°'^ oc mf,) and thus the gravity available to 
drive the collapse of the baryons prior to recombination on 
small scales like the third acoustic peak. This is highlighted 
in Fig [2] where the comparison is made between one sterile 
neutrino and two. 

The parameters for the best fit are given in table 1 
which also contains t he parameters for the WMAP5 fit from 
iDunklev et al.l (|2008l ) and a comparison of the two fits are 
shown in Fig [1] All parameters are consistent with exper- 
imental bounds and are not significantly different to the 
ACDM model, which is sensible since the ACDM model 
of the CDM anisotropies is a good one. 

The mass of the sterile neutrinos infered from the 
best fit value of Q^^h'^ = 0.117 is ~ lleV. This 

mass range of sterile neutrino has never before been con- 
sidered in the literature because it is excluded by cos- 
mological data if we assume Newton's law are correct 
jDodelson et al.l l2006l : ISeliak et all \200dt ) since they can- 
not influence galaxy rotation curves becaus e they would 
have a free streaming scale (cf. ISandersll2007h of more than 

= 1-3 (t?f)"*^' im^Y^' = SOfepc in a Milky Way 
type galaxy, for Vr — 200 km s ^ . The total mass this would 
create within 8kpc is ~ 5 x IO^Mq which is about 10% of the 
total mass and would actually help MOND fits to the Milky 
Way's rotation curve l|Famaev fc Binnevll2005l ; lGentile et al.l 
l2008l : lMcGaughll2008l '). 
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Figure 1. Shows the data of the CMB as measured by 
the WMAP satelli te year five data release (filled circles, 
Dunklev et"al] |2008|) and the ACBAR 2008 llReichardt et alj 
2008h data release (triangles). The lines are the ACDM max like- 
lihood (dashed) and the solid line is the fit with an lleV sterile 
neutrino with paramters given in table 1. The for the ACDM 
model has been scaled from the quoted 0.963 in IDunklev et all 
l|2008l ) to 0.979 here to better match the data. 



In particular, it would have a similar contribution to 
the energy density as required from cold dark matter fits 
to the CMB (r2^>^ = 0.117; Q.cdmh^ = 0.108) and leave 
the matter power spectrum at large scales (> 50/i~^Mpc) 
unaltered. This is shown in Fig [3] which compares the ob- 
served matter power spectrum with that predicted by the 
sterile neutrino model here, but with Newtonian instead of 
MONDian gravity. At scales smaller than ~ 50/i~^Mpc the 
computed power spectrum drops many orders of magnitude 
below the observed one. 

Qualitatively, this discrepancy is owed to the fact that 
structures on these small scales have formed with the as- 
sistance of MONDian gravity. For instance, following the 
argument of Eq[T] the redshift by which scales as large as 
50/i~^Mpc are deep in the MOND regime (i.e. g ^ ^) is 
roughly 



2<7 



1/3 



(4) 



which for 70Mpc is z » 1 00. Certainly many a uthors 
(ISanderd l2008l : iNusserl I2OO2I : iKnebe fc GibsonI |2004|) have 
shown that structures can form very quickly in MOND even 
without CDM and galaxy size objects can be in place as 
early as z ^ 10. 

The tools to perform the full matter power spectrum 
analysis are currently not available for MOND (nor stan- 
dard dynamics), since they crucially depend on hydrody- 
namics. Assuming that including the modified dynamics en- 
ables a match to the matter power spectrum at all scales, the 
only conceivable ways of distinguishing between MOND and 
ACDM (if missing satellites, the lack of cusps in DM halos 
and tidal dwarf galaxies are ignored) is in the complex mod- 
elling of galaxy formation, or the unambiguous detection of 
the hot or cold DM particles. 
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Figure 2. As for Fig[T]with the solid line again the single sterile 
neutrino fit with parameters given in table 1, but the dashed line 
is the fit with 2 sterile neutrinos sharing and the dotted line is 
with the maximum active neutrino contribution with f2y = 0.136 
and Ha compensated for a fiat universe. We reduced to 0.856 
and 0.939 respectively to match the amplitude of the first acoustic 
peak, but the second and third peaks are badly matched because 
there is not enough neutrino DM density on small scales because 
of the Tremaine-Gunn limit. 
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Figure 3. The filled circles are the data points from the SDSS 
llTegmark et al.ll2004) . the solid line is the single sterile neutrino 
model but with Newtonian instead of MONDian gravity. The 
dashed line is the ACDM model. 



3 DISCUSSION 

A single massive sterile neutrino appears consistent with 
the current level of precision in the measurements of the 
CMB anisotropies. It is also consistent with the mat- 
ter power spectrum at large scales (> 50/i~^Mpc) and 
is able to clump together with densities surpassing the 
maximum density of the DM in groups and clusters of 
galaxies where MOND requir e s dark matter of some form. 
As discussed in lAngus et al.l (|2008l ) there appears to be 



Parameter 


Single Us 


ACDM 


Active 


Two Us 


Ho 


71.5 


72.4 


71.5 


71.5 


lOO^t/i^ 


2.4 


2.27 


2.4 


2.4 




0.117 


0.0 


0.695 


0.117 


O h'^ 


0.0 


0.108 


0.0 


0.0 




0.965 


0.979 


0.856 


0.939 


No. massless u 


3 


3 





3 


No. massive u 


1 





3 


2 



Table 1. List of parameters used in the figures. The ACDM 
numbers come from Dunkley et al. (2008) but ns has been scaled 
from the quoted 0.963 to 0.979 for a better match to the data. 



a scale at which MOND begins to poorly describe the 
dy namics of astrophy s ical s y stems. This is highlighte d 
by Romanowsky ct al. (2003); Mikrom fc SandersI lj2003l l: 
[Angus et al.„ (.2008); O'SuUivan et al.. (,200"7l ) which show 
that no dark matter is necessary to explain the detailed dy- 
namics of relatively low mass groups of galaxies and sys- 
tems smaller. This is expected for sterile neutrino dark 
matter because it would have a free streaming length 
greater significantly larger than a typical galaxy (~50kpc 
for the Milky Way. However, just as numerical simulations 
of clusters of cold dark matter were necessary to show 
that the CDM halos are a poor match to observed galax- 
ies Jde Blok fc McGaugblll998l: iMcGaugh fc de Blo5 ll998l: 



Gnedin fc Zhad I2OO2I : iGentile et al.l l2004l : ICilmore et al l 
20071) ■ the equilibrium distribution of the sterile neutrino 



DM must be checke d to be consistent with groups and clus- 
ters of galaxies (see |Sanders 2007). 

On the other hand, the three active neutrinos should 
probably have masses well below Q.5eV. Otherwise it will 
become difficult to match the CMB power spectrum because 
the angular scale of the peaks prefers = 0.117 while 

ill/ oc m^. Increasing the mass of another neutrino reduces 
the mass of the sterile neutrino and the amplitude of the 
third peak of the CMB diminishes due to the rapidly de- 
creasing maximum density (p™""^ tx mf.). 

Certain analyses of neutrino mixing experiments seem 
to require an additional, sterile neutrino with a mass in the 
range AeV < m^^ < 18eV. Here I took the ansatz that there 
is a fourth, sterile neutrino of lleV mass and that MOND 
effects are not important at cosmological scales. I showed 
that its contribution to the dynamics of galaxies would be 
negligible, but that it could solve all problems MOND has 
with the dynamics of clusters of galaxies and it can match 
the angular power spectrum of the CMB. The matter power 
spectrum needs to be recalculated because MOND gravity is 
crucial to the formation of these smaller structures and be- 
cause of the increased dominance of baryons at these scales 
over DM, hydrodynamics cannot be avoided as in CDM sim- 
ulations. If experiments can indeed pinpoint the existence of 
a sterile neutrino with mass ~lleV this would be a signifi- 
cant advance for the Modified Newtonian Dynamics. 

Even if collider experiments detect a CDM candidate 
with mass of 300GeV, this will give us virtually no informa- 
tion about the cosmological abundance and therefore brings 
us no closer to solving the dark matter problem. The great 
thing about sterile neutrinos is that if we can find the mass 
from laboratory experiments then this effectively fixes the 
cosmological abundance AND the contribution the neutrinos 
can make to clusters of galaxies can be strictly constrained. 
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Henceforth, it would be possible to run structure formation 
simulations in MOND with all the ingredients. 
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